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Abstract. This paper presents an analysis of lensing of 
unresolved background stars. Previous calculations of the 
lensing rates and optical depths considered only resolved 
stars. However, if a faint unresolved star lens is close 
enough to a resolved star, the event will be seen by the 
microlensing experiments and attributted to the bighter 
star. The blending biases the duration, making the con- 
tribution of the unresolved stars very significant for short 
events. It is especially annoying, because this contribution 
is confused with lensing by brown dwarfs. The exact rates 
of these blended events are extremly sensitive to the lim- 
iting magnitude achieved in the microlensing search. Ap- 
propriate calculations of the optical depth and rates are 
provided here, and illustrated in the case of the DUO and 
OGLE experiments. The additional contribution of unre- 
solved stars is very significant. It probably explains the 
high optical depth and rates observed towards the Galac- 
tic Bulge. The blended unresolved event can be identified 
using either the color shift or the light curve shape. How- 
ever, neither of these two methods is apropriate to iden- 
tify a large number of blended events towards the Bulge. 
In some case of good photometry and small impact pa- 
rameter an indentification is possible. This is illustrated 
by the case of the OGLE 5 event, which clearly appears as 
a case of lensing of an unresolved star. The recent results 
obtained by the PLANET collaboration indicate that a 
high resolution and dense sampling of the light curve is 
possible, and will probably provide a very interesting pos- 
sibility to correct the blending bias, as demonstrated for 
OGLE 5. This possiblity, is certainly better than a sta- 
tistical estimation of the lensing rates, which are always 
prone to some uncertainty. But, at this time, the analysis 
of Microlensing events found in the various microlensing 
experiments requires the uses of modelisations of the con- 
tribution of unresolved stars. 
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1. Introduction 

In previous microlensing models the rates and opti- 
cal depth per star are usually estimated (Kiraga and 
Paczyhski, 1994 Han and Gould, 1995, MoUerach and 
Roulet , 1996). Below some limiting magnitude these stars 
can not be resolved individually and form a dense stellar 
background on all the surface of the images. I will call 
these star "background stars". These background stars 
can be lensed, and this issue has already been investi- 
gated in the case of the LMC (Bouquet, 1993), and for 
the Galaxy (Nemiroff, 1994). However the stellar fields of 
the current Microlensing experiments MACHO (Alcock et 
al. 1996, Alcock et al. 1993), OGLE (Udalski et al. 1994), 
DUO (Alard 1995a, Alard et al. 1995b) , EROS (Aubourg 
et al. 1993.) are extremly crowded, so tha inside the res- 
olution radius of a monitored stars, there are plenty of 
these unresolved background stars. Thus if one of these 
background stars is lensed, we will see a blended event. 
The amplification, duration and impact parameter will be 
biased by the blending. We may now wonder if their con- 
tribution to the total number of events observed is sig- 
nificant, compared to those expected from the resolved 
stars. As the blending of a background star with a re- 
solved star will be high we expect a large reduction of the 
effective Einsten Radius (Di Stefano and Esin, 1995). Con- 
sequently the observed number of events will be a compe- 
tition between the decrease of the Einsten radius, and the 
increasing number of stars at larger magnitudes. Also the 
effective duration of the event will be much shorter than 
it is for to events on resolved stars, and consequently it 
may produce a tail of short events. These short events 
could be confused with low mass lens events, and espe- 
cially lensing by brown dwarf, so that a detailed study 
of these events is particulary important. Another problem 
will be to recognise these blended events. The blending 
of light produces a modified light curve which might be 
distinguished from the unblended light curve. However to 
make an unambiguous separation a good time sampling is 
required, and a very accurate photometry is required. In 
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the case of these short events, we expect that it will be 
difficult to achieve a sampling dense enough. 



2. The Bulge luminosity function 

The study of lensing of unresolved stars is very sensitive 
to the shape of the adopted luminosity function. In previ- 
ous works (Kiraga and Pazcynski. Roulct and MoUcrach. 
Han and Gould), the number of stars at a given luminos- 
ity L is expressed as: dN = L~", with usually a ^ 2. A 
more realistic luminosity function can be derived from the 
Holtzman (Holtzman et al. 1993) observations of a field in 
Baade Window with the HST. The final Holtzman et al. 
luminosity function can be approximated by three con- 
tinuous straight segments with different slopes. The first 
change in the slope occur at F ~ 20 . With the extinction 
and distance to the Galactic Center adopted by Holtzman 
et al. , it correponds to a value of My = 4.9. . The Holtz- 
man et al. study agrees well with Terndrup (Terndrup 
1989) investigations of the Galactic Bulge which found a 
value of Mv = 4.07 for the Bulge turn-off at 6 = -8. In 
the Holtzman luminosity function the turn-off is situated 
around My — 4.0 also. The second change in the slope 
occur around V = 21.5 and seeems to continue until the 
Holtzman et al. limit situated at F ~ 22.5. A compar- 
ison with the Luminosity function of the globular clus- 
ters reviewed by Mould (Mould 1996) indicates that the 
shape of their luminosity function is almost the same in 
this range of magnitude. The Mould diagram allows also 
to a better determination of the last segment of the lu- 
minosity function. This straight segment seems to hold 
until My = 11 where a possible turnover in the lumi- 
nosity function is observed. However My = 11 is about 
7 magnitudes fainter than the limiting magnitude of the 
current microlensing surveys, and such faint unresolved 
stars do not contribute to the optical depth. The duration 
of the events on stars blended with a star 7 magnitudes 
brighter would be so short that it is not observable. Con- 
sequently the exact behavior of the luminosity function 
in this region is not important for our study, and we will 
make the reasonable choice to cut the luminosity function 
at My = 11. Looking now at the bright side of the lumi- 
nosity function, we find that at magnitudes brighter than 
the turnoff , the slope is very steep indicating that the 
number of stars drops rapidly. In this region the Bulge 
stars are essentially sub giants. The density rises again in 
the clump giant region, however these stars do not makes 
more than 5 percent of the stars in the current surveys. 
Consequently they have not a very significant influence 
on the total luminosity function. These stars should be 
treated separately (Alard 1996), and in the present study, 
we will ignore them. This will put the bright end of the 
luminosity function at My ~ 3 in the Holtzman et al. di- 
agram. The adopted value for the distance to the Galactic 
center in this study is 8.5 Kpc. The final adopted lumi- 
nosity function is shown in figure 1. 
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Fig. 1. The adopted Luminosity function for the present study. 

3. Estimation of the lensing rates and optical 
depth 

It is clear that, in an investigation of the lensing of unre- 
solved stars, consistency requires that we calculate the 
rates and optical depth also with the same luminosity 
function. This involves some slight changes from previous 
analysis. Therefore I will first investigate the case of re- 
solved stars considering a more general form for the lumi- 
nosity function. I first calculate the change in the optical 
depth. 

3.1. Optical depth 

The Numbers of source stars at distance Dg, and with 
absolute luminosity L can be written as: 

dN{Ds,L) = k lf{L) n{Ds) dL dDs 

Where k is a constant. 

Suppose that the experiment is able to find a total of Ntot 
stars up to a limit in apparent luminosity Iq. This leads 
to the following expressions: 

^^^l^IlB. = 1 X IfiD nm Dl dL dD, (1) 

with: 

Id= / lf{L) n{Ds) Dl dL dDs 

and: Lmi„ = IqX Dq/D^, where Iq is the apparent limiting 
luminosity. 

Lmax = hat X D^/Dq. Where hat is defined as the sat- 
uration limit of the detector. It can be expressed as: 
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hat = d X lo, where d is the detector dynamic range, a 
reasonable vahic is about 6 magnitudes in crowded fields, 
either for CCD or photographic plates. 

The optical depth associated with these stars can be 
calculated as (Paczynski 1991, Kiraga and Paczyhski 1994, 
hereafter refered as KP): 



droiD.,L)='J^^.fiD.) 



(2) 



with: 



f{Ds) = 



4ttG 



dDd 



£)2 

and I — L X 

To get the total optical depth we have now to integrate 
eq (1) over the source distances Dg and the absolute lu- 
minosity L of these stars. 

TO = ^ / Nlo{Ds) f{D,) n{D,) D^-^^ dD, 

-'fc JD, 



with: 



and: 



Nlo{Ds) = -^x Df / lf{L) dL 



L=L„ 



n{Dg) I?^2-2/3 



Here the factor is put in this expression to allow di- 
rect comparison with the standard optical depth formulae 
(Kiraga and Paczyhski 1994): 



) n{D,) D,^-^^ dD, 



We see that the Kiraga and Paczyhski expression is sim- 
ply cquivalant to taking Nlq{Ds) = 1 at all distances. 
A slight improvement is to take Nlo{Ds) = beyond a 
given distance (Roulet and Mollerach, 1995). The plot in 
Figure 2 allows a direct comparison between Nlo{Ds), 
and the approximation Nlo{Ds) = 1- Figure 2 illustrates 
the calculations both for DUO and for the OGLE experi- 
ment, which is about one magnitude deeper. Let's see now 
the changes in the total optical depth introduced by the 
Nlo{Ds) function. The optical depths are computed using 
the COBE bar model for the Bulge (Dwck et al, ..), and 
the Bahcall model for the Disk (Bahcall and Soneira, 1984) 
The ratio of the modified optical depth to the standard 
one is shown in table 1. The change is already important 
for OGLE, and is more significant again for DUO. It is 
easy to understand that the lower optical depth of the 
new formulae is essentially due to the fact that due to the 
limiting magnitude, the experiments are more sensitive to 
sources in the near end of the bar, for which the density 
of lenses on the line of sight is smaller. 



5 10 15 20 

Distance (Kpc) 

Fig. 2. Plot of the Nlo{Ds) function as a function of distance. 
For the DUO experiment a limiting magnitude of My = 3.5 

at the distance of the Galactic Center was adopted, OGLE is 
assumed to be one magnitude deeper. The dotted line is for 
DUO, and the dashed line for OGLE. 



Table 1. Optical depth ratio for DUO and OGLE 



Optical depth ratio 


Bulge lenses 


Disk lenses 


Total 


OGLE 


0.85 


0.89 


0.86 


DUO 


0.75 


0.8 


0.77 



3.2. Event rates. 

The Formulae given for the lensing rates by Kiraga and 
Paczyhski is: 

Tfe = / MDs) n{Ds) D/-^^ dDs 



Where: 



f^{D,) = -^l I p^{Dd)V{D,,Dd) 



DdjDs - Dd) 



1/2 



dDd dVi 



V{Ds, Dd) is the relative velocity between lens and source, 
and the Vi are the four components of this velocity. 
Now, it is straightforward to introduce the Nlq{Ds) func- 
tion in this expression, exactly as was done for the optical 
depth. 

This leads to: 

Fo = ^ / Nlo{Ds) MDs) n{Ds) D^-^^ dD, 
-'fe Jd, 
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The result of the calculation of the Fq and Tk rates is 
shown in figure 3 for DUO, and in figure 4 for OGLE. 
For all the figures the velocity dispersion given by Han 
and Gould (Han and Gould 1995) was adopted. The mass 
function is a Salpetcr mass function with a lower cut off 
at O.OSMq and an upper cut off at IMq. 
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Fig. 3. Comparison of the total rates of events for Bulge 
sources using KP formulae (dashed line), and the formulae de- 
scribed in this text (solid line). The total rate is estimated for 
the 1994 season of the DUO experiment using the DUO effi- 
ciencies (Alard and Guibert 1996). The two bumps are due to 
the particular shape of the DUO efficiencies. 



4. Modelling the lensing of blended unresolved 
stcirs 

4.1. Basic Principles 

The Einsten radius is defined as the distance for which the 

amplification of a microlcnsing event is a factor of 1.34. 
In the case of an unresolved background star of luminos- 
ity blended with a resolved star of magnitude Lq, the 
background star amplification required to make a total 
amplification in the combined image of 1.34 is: 

Ab = 0.34 X fb + 1.34 

Hereafter, I will refer to /& as the blending factor: 

fb = Lo/lo 

This gives the following expression for the blended Einsten 
radius: 

Rbifb) = ReX rcbifb) 




J I I I \ L 

20 30 
duration 



40 



50 



Fig. 4. Comparison of the total rates of events for Bulge 
sources using KP formulae (dashed line) , and the formulae de- 
scribed in this text (solid line). The total rate is estimated for 
10® stars in Baade window using the OGLE efficiency (Udalsky 
et al. 1994) 



with: 



rebifb) 



2{-l + Abl^{Al-l.Q) 



1/2 



4-2. Optical Depth due to the Unresolved stars. 

Let us consider unresolved background stars with absolute 
luminosity L at a distance Ds, whose number is dN(_Ds,L), 
and let us suppose also that they are blended with stars of 
apparent luminosity lb. The optical depth per resolved star 
associated with these stars can be calculated using equa- 
tion (2) by replacing the Einsten radius by the blended 
Einsten radius. 



dT{Ds,L) 



dN{Ds,L) 
N,„, 



X rb{lb,l) X f{Ds) X R 



The expression for dr contains the factor R, which takes 
into account that only those unresolved objects which are 
lensed close to a resolved star will be seen. The Value of 
R will be close to the ratio of the area of the seeing disk, 
to the mean space occupied by a star on the image, 
with: 

''6(^6,0 = [rebifb)?, and fb = h/l- 

We have now to take into account that the luminosity 
lb of the resolved star is not constant, but comes from the 
luminosity function. It was shown by Zhao (Zhao 1995) 
using the OGLE data that the number of resolved stars 
in Baade'swindow can be well described by a power law, 
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with an exponent close to —2. To be precise we have a 
fraction of stars near luminosity h: 



dF{h) = p2 lb"' C dh 



(3) 



with: a2 — 2 and /32 
consequently, 



Q!2 



dT{Ds,L,h) = dT{Ds,L) X dF{h) 

To get the total optical depth associated with the unre- 
solved stars, wc have now to integrate over the luminosity 
range of these stars, from the absolute luminosity (Lmm) 
to the faint end of their luminosity function (Lend)- 

In addition, we have also to integrate over all the possi- 
ble distances for the sources, and all apparent luminosities 
for the blends. 



This leads to: 



dN{D„L) 
Ntot 



xdF{lb) R n{lb,l) f{Ds) dD, dL dh 
If now we uses eqs, (1), (2), (3), we can express r as: 



S2 



20 40 

Distance (Kpc) 

Fig. 5. Plot of the Nl{Ds)/R function as a function of dis- 
tance for unresolved sources. The limiting magnitude defined 
for figure 1 were kept for this illustration. The dotted line is 
for DUO, and the dashed line for OGLE. 



R 



T = 



X lf{L) f{D,) n{Ds) Dl dL dh dD, 

This expression can be integrated over the variables L and 
4, consqucntly: 

T = ^ / Nl{D,) f{Ds) n{D,) D^-^^ dD, 



with: 



xlf{L) dL dh 

Note that the dependance in distance is hidden in the vari- 
able I = L X Dq/D^, and also in the boundary Lmin- 

Figure 5 illustrates the new Nl{Ds) function for the un- 
resolved stars in the case of the DUO and OGLE exper- 
iments. Note that these functiond reach their maxima at 
larger distances than the previous Nlo{Ds) functions, it 
means that a large contribution will come from unresolved 
sources on the far side of the bar. It is now possible to 
quantify the contribution of the unresolved sources to the 
total optical depth. In table 2 the ratio of the optical depth 
of unresolved sources to the the optical depth of resolved 
sources is indicated for DUO and OGLE. The calculation 
of this ratio require an estimation of the R constant. A 



{h,l) 



crude estimate of R is the ratio of the surface covered by 
the resolution radius to the mean area occupied by a star. 
For the DUO experiment the resolution radius is close to 
3 pixels, and the mean area occupied by a star is about 60 
pixels (Alard and Guibert 1996). This gives R ~ 0.5. For 
OGLE the resolution radius is probably close to the see- 
ing value (this is due to smaller pixels). Consequently with 
a mean seeing of 1.25" in Las Campanas, and a pixel of 
0.44" we can estimate the resolution radius as 2.7 pixels. 
The OGLE experiment follows 1.3 10^ stars on 14 images 
of 2048 X 2048 pixels each. This gives again a value for R, 
of ~ 0.5. 



Table 2. Ratio of the optical depth of unresolved sources to 
the the optical depth of resolved sources . A value R = 0.5 
is adopted for the calculations (see text for discussions). 



r/ro 


Bulge lenses 


Disk lenses 


Total 


OGLE 


0.58 


0.55 


0.57 


DUO 


1.93 


1.81 


1.89 



6 



C. Alard: Lensing of unresolved stars towards the Galactic Bulge. 



4-3. The rates of microlensing events from unresolved 
stars. 

The rates per resolved star for unresolved events can be 
calculated in the same way as for the Optical Depth. 
We have: 

r = ^ /" Nlj{D,) f^{D,) n{D,) Dl-^f' dD, (4) 
-'fe JDs 

with: 

xlf{L) dL dlb 

It is clear that the change of the Einsten radius by the 
factor reb{li,, I) changes also the duration by the same fac- 
tor. Consequently in the differential rates calculation, each 
duration ts = rs/V will have to be scaled by the factor 
rej,{li,,l). However, we see in equation 4 that for a given 
Einsten radius there is a distribution of the scaling factors 
reb{lb, I) hidden in the double integral -/Vi-y(_D.,). This dis- 
tribution can be calculated for each distance Dg, and for 
the realisation of the calculations, the scaling factor of the 
Einsten radius will be choosen by a Monte-Carlo method 
using this distribution. Some examples of such distribu- 
tions for different distances are shown in figure 6, in the 
case of the DUO experiment. 

The result of the calculation of the Fq and rates is 




0.0 0.2 0.4 0.6 0.8 1.0 

Duration scaling factor 

Fig. 6. Examples of duration scaling factor distribution func- 
tions, for different distances. The solid line is for a distance of 
8 Kpc, the long dashed line is for 15 Kpc, and the short dashed 
line is for 6 Kpc. 

shown in figure 7 for DUO, and in figure 8 for OGLE. It 
is evident that the DUO experiment is dominated by un- 
resolved stars. This result is in very good agreement with 
the data (Alard and Guibert 1996). Even for OGLE which 



has the deepest photometry, the bias is still very signifi- 
cant. It is rather straightforward to relate these events to 
the excess in the rates and optical depth observed towards 
the Bulge. This excess was not been explained, even us- 
ing a bar in the Galactic model. We find here a natural 
explaination. 








10 20 30 



duration 

Fig. 7. Comparison of the the contribution to Bulge microlens- 
ing event rates for resolved stars (short dashed line), and unre- 
solved stars (solid line). The total rate is represented by a long 
dashed line. The rates are estimated for the 1994 season of the 
DUO experiment using the DUO efficiencies (Alard and Guib- 
ert 1996). Note the dominant contribution of the unresolved 
stars. 



5. Can we identify the blending events using the 
light curves ? 

The light curve of a blended event is affected by the addi- 
tional light coming from the companion, so that the light 
curve is modified. This blended light curve has more ex- 
tended wings than an unblended event, so that it might be 
possible to differentiate it from the unblended one using a 
chi-square test. 

Another important issue is the possibilty that the ampli- 
fied star and its companion may have different colors. For 
instance if the companion is more blue than the ampli- 
fied star, we expect more blending in the blue, and thus a 
larger reduction of the amplitude in this color than in the 
red. 

5.1. Color variations 

Let us now investigate the problem of color variation dur- 
ing the event in more detail. In the current microlensen- 
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Fig. 8. Comparison of the contribution to Bulge microlensing 
event rates for resolved stars (short dashed line), and unre- 
solved stars (solid lino) . The total rate is represented by a long 
dashed line. The rates are estimated for 10*^ stars in Baade 
window using the OGLE efficiency (Udalsky et al. 1994). 



ing experiments the sources are concentated in the range 
18 < V < 20. An examination of Tcrndnip (Tcrndnip, 
1988) color magnitude diagram in Baade'swindow indi- 
cates that the mean V — I color is about 1.3 to 1.4. For 
fainter stars, where expect to find the unresolved sources, 
the Holtzman et al. color magnitude diagram shows that 
the color changes again very slightly. The mean color a 
close to 1.4 at V = 22. We sec that in the apparent mag- 
nitude range of interest we expect a diferential color vari- 
ation of about 0.1 magnitude at most between the un- 
resolved source and the blending companion, we have to 
reach = 23 to expect more significant color variation 
between sources and blend (about 0.2 magnitudes). Un- 
fortunately the photographic technique does not perform 
very well in the red band, thus an accurate investigation of 
the color changes during the event is not possible. It means 
that the DUO survey will not be able to detect the slight 
color changes expected. The OGLE survey has a good cov- 
erage only in the I band, the sampling in V is very sparse, 
and does not allow color analysis during the event. In the 
MACHO case, photometric data are available in two large 
bandpasses, which is transformed to a,V — R color index. 
The analysis of the LMC events by MACHO shows that 
they can identify color variations due to blending for some 
of the events. However there is a small residual noise on 
the color variations of all the events of about 0.02 to 0.03 
magnitudes, associated with the photometric errors. At 
a 3 a confidence level we may require a color variation 
of about 0.075 magnitude to firmly established the chro- 



maticity of the event. Converted to V — 7 it requires a 

color changes of about 0.15 magnitudes at least to recog- 
nise a blended event with a good confidence level. 
We consider now with a such accuracy how many of 
the Bulge microlensing events which involves unresolved 
sourves could be identified as such. 

Let us make the following simple model: 

And we assume that the difference of color can be de- 
scribed by a gaussian distribution, shifted by a systematic 
value. Thus the number density n(C) of stars with color 
a difference in color C will be expressed as: 



n{C) = exp 



(C - shiftc)^ 
-2 a? 



(5) 



The Holtzman et al. color magnitude diagram indicates 
that for the gaussian we can assume ctc — 0.07. In the 
range U = 20 to U = 22, a conservative value for the shift 
is shiftc ~ 0.15. The color shift is then given by: 

Ac = 2.51og(^,/A,) 

i and v stand where the subscripts for the two photometric 
colors, with: 

1 + afb 



and: 



1 + fb 

1 + afbC 



1 + /6C 

a and c are defined by the following expressions: 



a = 



u^ + 2 



U X 



c = 10 ^ 



The number of events with color blending signature iden- 
tified will be simply the number of events with a color shift 
Ac > 0.15. It can be easily calculated if we assume that 
the unblended impact parameter u has a uniform distri- 
bution. In this case it is sufficient to integrate the number 
of events with Ac > 0.15 over the u distribution with a 
given color difference C. The calculation is completed by 
suming over the color difference distribution defined in eq 
(5). The maximum value of the impact parameter is set by 
the blending factor fb using the formulae: 



^max — 



with: 



2{-l + A,/^{Al-l.Q) 



1/2 



(6) 



Ab = 0.34 X h + 1.34 

The calculation of the percentage of events with color 
shifts is performed for several values of the blending factor 
/(,. An illustration of these calculations is given in Figure 
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9, where the color shift is expressed as a function of the im- 
pact parameter for a few fb values, and color diffrenee of 
0.2 magnitude. Looking at this diagram we can already 
guess that the number of events with detectable color 
shitfts will be small. The final values of the events with 
color shifts identified is given in table 3. The values are 
extremly small. A shift of 0.15 mags in the color difference 
distribution, even for high values of the blending factor. 
The values for a shift of 0.3 magnitudes are also given, 
to illustrate the case of the very faint unresolved stars. 
However these stars will be extremly blended, and con- 
sequently will give events will very short duration which 
will be almost all removed by the low efficiency of the 
experiments in this range. Consequently the color shift 
method will give poor results in the case of the Galac- 
tic Bulge. In the case of the Magellanic Clouds, the color 
changes rapidly with the magnitude close to the limit 
of the microlensing experiments, which explains why so 
many events with color shifts are found by MACHO (Al- 
cock, et al., 1996)towards the LMC. 




0.2 0.4 0.6 0.8 



(u/umax) 

Fig. 9. The distribution of the color shift for blending pa- 
rameters of 2,5,10,50, and a color difference of 0.2. The curve 
with fb = 50 is upper, and the values 10,5,2, are below in a 
decreasing order. 



Table 3. The fraction of events with detectable color shifts. 



Blending factor 


2 


5 


10 


50 


shiftc = 0.15 


0.0025 


0.0219 


0.0364 


0.0485 


shiftc = 0.3 


0.1677 


0.3079 


0.3559 


0.3906 



5.2. The shape of the light curves. 

I will investigate another possibility to identify a blended 
microlensing event based on the shape of the light curve 
in this section. The blending of light modifies the shape 
of the light curve, and consequently we may use a statis- 
tical test to see if a light curves differs significantly from 
an unblended one. However we have to realise that the 
likely difference due to blending is rather small, and we 
may expect that with the current available photometry 
it will be difficult to identify a blended light curve at a 
significant level of confidence. It is possible to adress this 
problem in more general terms using Monte-Carlo simula- 
tion of microlensing events. Assuming a noise distribution, 
it is easy to simulate microlensing events with different 
blending factors. These blended events can be analysed 
by fitting an unblended curve to the simulated data set. 
We expect a systematic difference in the chi-square of the 
fir due to the different shape of the blended light curve, 
especially in the wings. However the problem is to know 
how significant is this difference compared to the normal 
chi-square variations for unblended events due to noise. 
This problem can be easily addressed if we are able to 
build the chi-square distribution for a series of blended 
light curves. Such chi-square distributions are illustrated 
in Fig. 10, the line shows the limit within which 95 percent 
of the chi-square distribution for unblended events is con- 
tained. Beyond this limit we have a 95 percent confidence 
level that the chi-square indicates a light curve which is 
systematically different from the unblended model. It is 
now sufficient to sum the fraction of the blended distribu- 
tion beyond this limit to get the fraction of blended events 
Rj, which can be identified with a good confidence level. 
The ability to recognise a blended event will of course de- 
pend on the amplitude of the event, for a given blending 
factor fb and a given duration to. The amplitude is related 
to the impact parameter u. Consequently to get the ex- 
pected fraction Ri, , at fb and to it is sufficient to sum on 
the impact parameter u, in the range to Umax ( Umax is 
defined in eq 6). The result of the corresponding calcula- 
tion is illustrated in Fig 11, as a function of the blending 
factor, for a series of durations. The blended events are 
well identified only for small impact parameters. Beyond 
Umax/W, the efficiency drops significantly, which makes 
the total efficiency quite low. These tabulated expressions 
of Rb as a function of fb and to are now introduced directly 
in the calculation of the unresolved event rates. Then new 
rates corrected for the number of blended events which 
could be recognised are computed. Fig 12 shows a com- 
parison between this rate and the uncorrected rate in a 
case resembling to the OGLE experiment. Typical errors 
distributions are extrpolated from the DUO experiment 
(Alard and Guibert 1996), and are divided by a scaling 
factor of 2 to take into account the better quality of the 
OGLE photometry. The resulting errors distribution has 
a sigma about 0.07 magnitude for most of the points, and 
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Fig. 10. An example of simulated Chi-square distributions for 
unblended events (full line), and for blended events (dashed 
line). The vertical line indicates the 95 percent limit for the 
unblended distribution (see text for explanations) . For this ex- 
ample the events are simulated with a duration of 40 days, and 
an impact parameter of 0.05. The blending factor has a value 
of 4. 




Blending Factor 

Fig. 11. The fraction of blended events identified as a function 
of the blending factor for different durations. The durations are 
respectively 100,63,25,15,10,6,4,2.5 days (from top to bottom). 



reaches 0.1 or slightly more in 10 percent of the cases; it 

seems to be an acceptable description of the OGLE pho- 
tometry on stars close to the limiting magnitude, which 
represents most of the sample. We see that the difference 
is small, and considering that OGLE has the best photo- 
metric accuracy, it is certainly worse again for the others 
experiments. 

The conclusion is again, as for the color test, that given the 
mean quality of the photometry available in the microlens- 
ing experiments only, a slight percentage of blended events 
should be recognised on the basis of the light curve shape. 
Only the the events with good photometry and small im- 
pact parameter should be identified as blended. However 
for some of the monitored stars the OGLE experiment is 
able to perform photometry much better again than for 
the majority of the stars. This gives an interesting oppor- 
tunity to look for a blending signature in the light curve. 



0.15 - 
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Fig. 12. The event rate statiscally corrected for the fraction 

of blended which might be identified on the basis of their light 
curve shape. The corrected rate is represented by a full line, 
and the initial rate by a dashed line. 



6. OGLE 5: lensing of an unresolved stcir 

The photometry achieved on the OGLE 5 event is of an 
interesting quality, the errors bars are as good as a few per- 
cent on many points. However this event is not well fitted 
with the standard unblended model (see figure 13). The 
discrepancy is especially large in the wings of the event, 
so that it is natural to try to fit this event with blending. 
Fig. 14 shows the dramatic improvement of the chi-square 
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when a blended light curve is fitted. The best fit is ob- 
tained for a blending factor of 2.45 (Fig 15.). Thus makes 
the lensed source 1.34 magnitude fainter, and places it 
below the OGLE detection limit. Consequently it is very 
likely that OGLE 5 is an example of lensing on an unre- 
solved star. The better Chi-Square per degree of freedom is 
different from unity, but we have to recall that an event on 
an unresolved star is seen only by blending on a resolved 
star. For the photometry, the position of the resolved star 
is used, but it might be significantly different from the 
position of the true magnified star. The achieved errors in 
a PSF fitting routine are rather sensitive to the quality 
of the positioning, thus we expect somewhat larger errors 
in the case of an unresolved event. The small discrepancy 
in the Chi-Square per degree of freedom is thus perfectly 
consistent with the scenario of an unresolved event. 
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Fig. 13. The fit of an unblended light curve to the data, note 
the large discrepancy in the wings. 



7. Discussion 

In this paper I show that the rates of microlcnsing events 
and optical depth to microlcnsing due to the unresolved 
stars can be modelled. However I emphasize that the ex- 
act rates and optical depth arc closely related to the num- 
ber of unresolved stars which will contribute per resolved 
star. This number is well constrained in crowded fields, 
essentially because the surface occupied by a star on the 
image is set by the crowding limit. On the other hand, 
this crowding limit is closely related to the resolution ra- 
dius, thus the surface occupied by a star is just a function 
of this resolution radius. A value of 0.5 was found for R 
(the ratio of the surface covered by the resolution radius 
to the mean surface occupied by a star), both for DUO 
and OGLE. This simply means that probably there is a 
linear relation between the resolution radius and the mean 
radius per star. 



Fig. 14. The Chi-Square per degree of freedom obtained for a 
fit to the data for OGLE 5 with different blending factors. 
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Fig. 15. The best fit with blending, note the large change in 
the estimated duration in comparison with the unblended fit. 



However in the crowded fields, if a new star appear at 
more than a resolution radius, it does not mean necessar- 
ily that it will be resolved. There are so many stars that 
it might well be confused with another one. This issue 
should be clarified later using Monte-Carlo simulation of 
crowded fields. This effect will probably lead to a slight in- 
crease of R. Another problem is that if the amplified star is 
situated at some distance from the resolved star the pho- 
tometry which assumes fixed positions (except for DUO) 
is certainly affected. This possibilty was already discussed 
in the case of OGLE 5 and will cause a slight drop in the 
efficiency, which will increase as the unresolved star will 
be more distant form the resolved one. This issue could 
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be also clarified later using Monte-Carlo simulation. But 

globally, this effect will more or less canccll out with the 
confusion effect exposed just before. Thus to conclude, the 
value R = 0.5 might be an acceptable approximation. 
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Fig. 16. Comparison of the total microlcnsing rates (unrc- 
solvcd+rcsolvcd stars) with the rates for resolved stars for dif- 
ferent cut-offs in the mass function in case of the OGLE exper- 
iment. The cut-offs are: 0.02 MQ (long and short dashed line), 
0.04 MQ (long dashed line), 0.06 MQ (dotted line), and 0.08 
MQ (short dashed line). The scaling factors are respectively: 
1.09, 1.24, 1.18, 1.28 (see text for explainations) . 



case of the DUO experiment will be treated in another pa- 
per (Alard and Guibcrt 1996). While the unresolved star 
bias has to be taken into account, as demonstrated in the 
previous section, the modelling includes some uncertain- 
ties. However for the data set already assembled by the 
microlensing experiments, it is certainly an acceptable so- 
lution to the unresolved star bias. For the future, a very 
interesting possibilty is offered by the PLANET (Sackett, 
1995) collaboration, who will provide a dense and accu- 
rate sampling of the events. It will allow a much better 
control of the bias caused by the unresolved, using both 
color shift but importantly the shape of the light curve, 
as demonstrated for OGLE 5. 

To conclude, I will also emphasize that the contribution 
of the unresolved stars drops if the magnitude limit is in- 
creased. This is well illustrated by the comparison OGLE 
vs. DUO. Consequently a first solution is to try to reach 
stars as faint as possible to minimise the bias. It is cer- 
tainly a possible orientation for the OGLEII experiment, 
which will achieve improved resolution better again with 
a new telescope. On the other hand for MACHO and 
EROSII experiments, which use bigger pixels, and cover 
larger fields, an interesting solution is certainly to moni- 
tor the bright clump giants. These stars are much brighter 
(about 3 magnitudes) than most of the bulge stars, and 
their density per image is also much lower. Thus, the 
imrcsolvcd stars per giant will be low and they will be 
so blended that the effective Einsten radius will be very 
small. This idea leads to the conclusion that the unre- 
solved star bias will be very small on clump giants. Some 
other advantages of the clump giants include the good 
photometry expected, and also the possibility to have a 
sample with a very high completeness (Gould 1995). 



8. Conclusion. 

The modeling of the Unresolved star microlcnsing rates 
demonstrates that they induce an important bias for the 
short events. I show that given the current photometric 
errors, most of these events would be hard to distinguish 
from the unblended events. Thus this is an important and 
annoying bias. If this bias is not taken into account, it 
will certainly influence the shape of the lens mass func- 
tion estimated from the data. For instance, if a Salpeter 
mass function is used, we expect that the lower cut-off 
will be shifted towards the brown dwarfs. This idea is 
illustrated in Fig. 16 where the total OGLE rates (re- 
solved-|- unresolved) are compared to the rates for Salpeter 
mass functions with different lower cut-offs (I recall that 
initially I used a Salpeter mass function with a lower cut- 
off of O.O8M0 to compute the OGLE and DUO rates). A 
better agreement with the total rates is found if the differ- 
ent trial models are slightly scaled. I calculated this scaling 
factor, so that the trial models give the same total rates in 
the range to 80 days as the total rate (see Fig 16.). The 
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